THE ZETA FUNCTION OF THE LAPLACIAN ON CERTAIN 

FRACTALS 
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Abstract. We prove that the zeta-function (a of the Laplacian A on a self-similar 
fractals with spectral decimation admits a meromorphic continuation to the whole complex 
plane. We characterise the poles, compute their residues, and give expressions for some 
special values of the zeta-function. Furthermore, we discuss the presence of oscillations in 
the eigenvalue counting function. 



1. Introduction 

Diffusion on fractals has been studied extensively as a generalisation of usual Brownian 
motion on manifolds. After its introduction in the physics literature (cf. |30j) M. Bar- 
low and E. Perkins [3] gave a very detailed study of the properties of the diffusion on 
the Sierpihski gasket. Later T. Lindstr0m generalised these results to nested frac- 
tals. Especially, he derived results on the distribution of the eigenvalues of the Laplacian 
associated to the diffusion. Here the Laplacian is seen as the infinitesimal generator of 
Brownian motion. Alternatively, the Laplacian can be obtained as the limit of properly 
rescaled difference operators on graphs approximating the fractal structure (cf. [3 13 HE] ) • 
Furthermore, the Laplace operator can be obtained via Dirichlet forms (cf. pQ). 

There exists a vast literature on properties of the spectrum of the Laplacian. For in- 
stance, we refer to UHl HH1 HZl HH1 EH I2H E21 EH EH1 EHI Especially, it has been proved 
that the Laplacian on an infinite post-critically-finite (p. c. f.) self-similar fractal has pure 
point spectrum (cf. [32J). 

In the case of a compact self-similar fractal G equipped with the Hausdorff measure 7i 
it is natural to ask for the behaviour of the counting function of the eigenvalues under 
Dirichlet or Neumann boundary conditions 

(1.1) L(x)= 

— Au=fiu 
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the trace of the heat kernel Pt(x, y) 

(1.2) P(t) = e ~ M = / Pt(x,x)dH(x), 

— Au=[MU 

or the corresponding Dirichlet generating function 
(1-3) Ca(«)= £ ^ 

the "zeta-f unction" of the Laplace operator (the eigenvalues are counted with multiplicities 
in all three sums). This zeta- function will be the subject of the present paper. 

In the case of the Laplacian on a Riemannian manifold, the meromorphic continuation 
of Ca can be derived directly from the asymptotic expansion of the trace of the heat-kernel 
(cf. fI7\ 13 lj). In this case it is known that admits a meromorphic continuation to the 
whole complex plane with at most simple poles at §, f — 1, . . . , I — L^"J > where d denotes 
the dimension of the manifold. Furthermore, the residues at these poles have a geometric 
meaning. 

In the case of a fractal the situation is more delicate. The current knowledge on the 
behaviour of the heat-kernel is not sufficient to provide enough information on its trace to 
get the meromorphic continuation of £a- In the present paper we will present a method for 
the continuation of (a which is based on an asymptotic study of the solution of Poincare's 
functional equation. In Section El we prove several results relating the behaviour of the 
solution to properties of the according Julia set. This approach for the analytic continuation 
of Ca is restricted to a subclass of self-similar fractals which exhibit spectral decimation 
(cf. jSD]) with a polynomial map (see Section |2J). On the other hand it allows to attach 
a zeta function to every polynomial which has a certain behaviour under iteration. Such 
zeta-functions have been introduced in [321; we give a totally different approach for their 
meromorphic continuation and are also able to derive expressions for some special values. 

In Section El we use the properties of the zeta function to confirm a conjecture (cf. j!7|. 
p. 105]) about the presence of oscillations in the eigenvalue counting function of the Lapla- 
cian for the class of fractals under consideration. 



2. Eigenvalues of the Laplacian and spectral decimation 

As it was observed in the case of the Sierpihski gasket and its higher dimensional ana- 
logues in EH EH] the eigenvalues of the Laplace operator on certain finitely ramified 
self-similar fractals exhibit the phenomenon of spectral decimation (cf. jU] for the case of 
higher-dimensional Sierpihski gaskets and [2H1EZ1I32] for a more general treatment). 

Definition 1 (Spectral decimation). The Laplace operator on a p. c. f. self-similar fractal 
G admits spectral decimation, if there exists a rational function R, a finite set A and a 
constant A > 1 such that all eigenvalues of A can be written in the form 

(2.1) A m lim \ n R(- n \{w}), weA,meN 
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where the preimages of w under n-fold iteration of R have to be chosen such that the limit 
exists. Furthermore, the multiplicities of the eigenvalues depend only on w and m, and 
the generating functions of the multiplicities are rational. 

The fact that all eigenvalues of A are negative real implies that the Julia set of R has 
to be contained in the negative real axis. We will exploit this fact later. 

In many cases such as the higher dimensional Sierpihski gaskets, the rational function 
R is a polynomial. The method for meromorphic continuation of (a given in Section El 
makes use of this assumption. In a recent paper ^2] A. Teplyaev showed under the same 
assumption that the (^-function of the Laplacian admits a meromorphic continuation to 
> — £ for some e > depending on properties of the Julia set of the polynomial 
given by spectral decimation. His method uses ideas similar to those used in ^U] for the 
meromorphic continuation of a Dirichlet series attached to a polynomial. Complementary 
to the ideas used here, Teplyaev's method carries over to rational functions R. 



3. Asymptotic behaviour of solutions of a functional equation 

In this section we collect several results on the growth of the solution of the Poincare func- 
tional equation related to polynomials. We note here that Valiron [33J Chapter VII] gave 
rather general results on the order of the solution of such functional equations. Throughout 
the paper we will use the notation 

logz 
logA 

for the logarithm to base A. We first give a coarse growth estimate, which is contained in 
[33] and will be refined in Theorem ^ 

Lemma 1. Let p(x) = adX d + • • • + a%x be a polynomial with real coefficients, d >2, which 
satisfies p(0) = and p'(0) = A > 1. Then the solution $ of the functional equation 

(3.1) $(\z)=p($(z)), $(0) = 0, $'(0) = 1 

is an entire function of order < p = log A d. 

Proof. By the general theory of iteration of polynomial functions, there exists a unique 
entire solution $ of the equation (|3.1j) (cf. jUEDJ)- Furthermore, there exists a constant 
C > 1, such that \p(z)\ < Cmax(l, \z\ d ). Therefore, 

|$(*)| < C cfl ~ 1+ - +(m max(i, m\- n z)\ dn ). 

Let r > such that |3>(;z)| < 1 for \z\ < r and take n = [log x (\z\/r)] + 1 to obtain 

\<$>(z)\ < C dr ~"^ P . 

□ 



The following theorem gives a more precise description of the growth order of $ under 
assumptions on the behaviour of p under iteration. 
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Theorem 1. Let p(x) = a^x d + ■ ■ • + a\x be a polynomial with real coefficients, d > 2, 
which satisfies p(0) = and p'(0) = A > 1. Furthermore, suppose that T^, the Fatou 
component of oo of p, contains an angular region of the form 

Wf, = {zeC\{0}\ |argz| <(3} 

for some (3 > 0. Then the solution $ of the functional equation (13. 1|) is an entire function. 
For any e > and any M > the asymptotic relation 

(3.2) $(z) = T -^exp (z»F(\og x z) + o (\z\- M )) 

holds uniformly for |arg^| < (3 — e, where p = log x d and F is a periodic holomorphic 
function of period 1 on the strip {w G C | \$sw\ < j^}- The real part of z p F(\og x z) is 
always positive; F takes real values on the real axis. 

Proof. After Lemma ^ it remains to prove only the assertion about the asymptotic be- 
haviour of 

Since p is a polynomial of degree d > 2, there exists a real number R > 0, such that 
\z\ > R =J> \p(z)\ > 2\z\, which implies that (p^ n \z)) n tends to infinity for \z\ > R. Fix 
e > 0. Then by local conformity of $ around there exists r > such that 

$ ({z G C \ {0} | |*| < r, | arg*| < (3 - e}) C Wp. 

By definition of the Fatou set the sequence (p^ n '(z)) n is equicontinuous on Wp and therefore 
contains a subsequence {p( nk \z))k, which converges to oo uniformly on 

M e = $ ({z G C \ {0} | j < \z\ < r, | arg z\ < (3 - e} 

by the fact that oo is an attracting fixed point of p and our assumptions on the Fatou 
component of oo. Thus there exists an N such that \p^- N \z)\ > R for all z G M £ , which 
implies that (p^ n \z)) n converges to oo uniformly on M £ . Inserting this fact into (j3.1|) this 
implies that 

(3.3) lim |$(-z)| = oo uniformly for | argz| < (3 — e. 

z— >oo 

Furthermore, for any K > 1 there exists a C > such that 

\®( z )\ > C\z\ K for | argz| < (3 — e. 

Now note that our assumption on the Fatou set implies that there are no zeros of $ in 
Wp. Therefore, 

ty(z) = log$(z) 
is analytic in Wp and satisfies the functional equation 



(3.4) -fy(Xz) = d^f(z) + log a d + r (z) with r(z) = log : , r 

^ rf ^ \Z) 

The function r(z) is analytic in 

{z G C | |*| > r , | &rgz\ < (3} 



k=0 
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for some r > 0. Setting ^f(z) = z p (j)(z) — with p = \og x d yields 

(3.5) <j>(\z) = fa) + 

By our previous knowledge on $(z) for z — > oo in | argz| < (3 it follows that r(z) tends to 
faster than any negative power of \z\ uniformly in any angular region | argz| < (3 — e. 
Iterating ()3.5j) yields 

(3.6) ^) = ^4E r ^ 
We set 

F(log A z)= hm <f)(X n z), 

n— too 

which exists by our knowledge on r(z) and ()3.6|) . Clearly, F(w) is periodic with period 1, 
holomorphic in \Qw\ < and 

<{>{z) = F (log A z)-^ r ~^0- = F (log, z) + o(\z\- M ) 

for any positive M. 

Putting everything together, we obtain 

9iz) = zf>FQ g x z)-'^ + o{\z\-»), 

which implies ()3.2)1 . The statement on the real part of z p F(\og x z) follows from |$(-z)| — > oo 
for \z\ — > oo in Wp. □ 

Corollary 2. The limit 

\og<&(\ n z) 
hm 

7i-*oo d n z p 

exists for all z G $^ 1 ^ ) (JF 00 ). This limit gives the analytic continuation of F(\og x z) to 

$ ( - 1} (^oo). 

Proof. The proof of the asymptotic formula for ^(z) used only the fact that $(A n z) — > oo 
for n — > oo. This relation holds for any point z G $^ _1 - ) (^ : ' 0O ). □ 

Remark 1. Constancy of the periodic function F occurring in the statement of Theo- 
rem[T]has been the subject of a series of papers in the context of branching processes 0IH], 
where it turns out that the constancy of F is implied by the existence of a continuous 
time extension of the according branching process (cf. ^3]). Furthermore, it is conjectured 
that also the opposite implication is true. Usually, on the real axis the function F exhibits 
very small oscillations around a mean value, which can only be observed using high pre- 
cision numerical computations. Therefore, theoretical conditions for the presence of such 
oscillations are of special interest. 
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For the purposes of this paper we give the following condition for non-constancy of the 
periodic function F. 

Theorem 3. Let $ be the solution of the functional equation (|3.1|) . Assume that there 
exists an angle 7 such that <£>({re 17 | r > 0}) intersects the Fatou- component Too as well 
as the Julia-set J v . Then the periodic function F in ()3.2)1 is not constant. 

Proof. We note first that J v is a compact subset of C, since p is a polynomial (cf. 0]). 
Assume that _F(log A z) = C. Then by Corollary |21 we have 

| . m log»(A !£ ) =c 

n-»oo d n Z p 

for z G $ (_1) (J' 00 ). For any r > with $(W 7 ) G J p we have $(rA™e i7 ) G J p for all n. 
Thus by our assumptions, for any fceN, any e > 0, and any M > there exists r > M, 
such that $(re 17 ) G .Foe and 



Then we have 



inf{|<l>(re i7 ) - z\ \ z e J p } <e. 



log$(A n+fc re i7 ) ^ 
nm ; tt. ■ = ^ 



l og p(fc)($(AVe i7 )) _ 
n-*x> d n (p( k \^(re i "')))P 



On the other hand 



lim lo g $ ( An+fere ' 7 ) _ jj logp^($(A n re i7 )) (p( fc )($(W 7 ))) p _ 

r p^n+fc e ip 7 " rfn(p(A:)($( re i7)))P d^Pe^P 



n— »oo 



Since p( fc )($(re J7 )) is contained in a bounded set and r can be chosen arbitrarily large. 
This gives a contradiction. □ 

Corollary 4. Assume that J p C M~ and i/iai $ zs unbounded on Mr. Then the periodic 
function F in Theorem^ is not constant. 

Proof. $(R~) intersects the Julia set J pi since $ attains negative values close to 0, and 
G J p . $(R _ ) intersects Foo, since $ takes arbitrarily large values. Then Theorem El gives 
the assertion with 7 = 71. □ 

Remark 2. The example of = x(4 + x), where <&(z) = 4sinh 2 \\fz-, shows that the 
assumption on the unboundedness of $ on R~ cannot be omitted. 

Corollary 5. If p < \, then the periodic function F in TheoremU\is not constant. 

Proof. By |UJ Theorem 3.1.5] a function of order < | is unbounded on any ray, and so is 
$. Apply Corollary 0] to a ray, whose <3>-image intersects the Julia set of p to obtain the 
assertion. □ 

Corollary 6. If the periodic function F in Theorem Q] is not constant, it has infinitely 
many non-zero Fourier-coefficients. 
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Proof. Since p(0) = and lies on the boundary of by the hypotheses of Theorem^ 
infinitely many zeros of $ lie on the boundary of (J 7 ^). From this we get that F(log A z) 
is unbounded when approaching these boundary points of $^ _1 ^(JF 00 ), which implies that 
F has to have infinitely many non-zero Fourier-coefficients. □ 

4. Spectral decimation and level sets of $ 

By (J2.1)) all eigenvalues of the Laplacian A can be computed by iterating p(~^ and 
rescaling. In this section we want to analyse this more precisely. The inverse map of the 
polynomial function p has d = degp branches, which we denote by qi, ... ,qd- We denote 
by qi the unique branch given by qi(0) = 0, which also satisfies q[(0) = v. Given w G C 
the elements of p(~ n '({w}) are in correspondence to words of length n over the alphabet 
{1, . . . , d}: associate to £ = (£i, . . . , £ n ) G {1, . . . , d} n the value z £ = qg n o q n _ x o • • • o q £l (w). 
Then p( n \ze) = w and all elements of p(~ n \{w}) are of this form. 

Let I = (l h £ 2 , . . .) G {1, . . . , 4 N Then the limit 

lim X n q ln o g n _! o • • • o q h (w) 

n— »oo 

exists if and only if there exists iV G N, such that £ n — 1 for n > N. This follows from the 
fact that q\ is the only branch of p^ 1 ^ which maps to itself. 

Take z n G X n p^ n \{w}). Then define u n G C by \~ n z n = $(A- n w„). Then z n = 
\ n $(\~ n u n ) and assuming that the limit lim^oo z n exists as discussed above, we have 

z = lim z n = lim X n §(\~ n u n ) = lim u n . 

n— >oo n— >oo n— >oo 

On the other hand w = p^(^(\~ n u n )) = $(«„) and therefore $(#) = w. Thus the 
solutions of $(2) = w are exactly the possible limits of sequences (z n ) as above. Summing 
up, we have 

Lemma 2. Let I = (l x , 4, . . .) G {1, . . . , d} n with £ N = £ N+1 = . . . = 1 for some JVgN. 
Then the limit 

Zl = lim \ n q ln o q n _ t o • • • o q h (w) 

ri —>oo 

exists and satisfies $(zg) = w. All solutions of the equation $(z) = w can be obtained as 
such a limit. 

5. Weierstrass products and Dirichlet series 

Since the solution $ of the functional equation (J3.1j) is of order p, we can represent 
®w( z ) = 1 — ^^(^) by an Hadamard product (for ro^O) 

(5.1) $ w (z) = exp( Cl (w)z + ■■■ + c k (w)z k ) 

(Z \ f Z Z^ \ 

_ 1 + M^yJ exp r^R + ^h 2 + ' ' ' + WM k ) ' 

where k = |_pj denotes the integer part of p (cf. jH]). It is clear from the validity of the 
functional equation ()3.1|) that <& w has infinitely many zeroes even in the case of integer p. 
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It follows from Theorem^ that log$ w (z) has a convergent Taylor series around 0: 



(5.2) log$ w (z) =J2blHz £ 



1=1 

Furthermore, assuming that the Julia set of p is a subset of the negative reals and w < 
as it is the case for spectral decimation, log<& w (z) is holomorphic in {z G C | | argz| < 0}. 
Computing the Taylor series of \og§ w (z) around z = and comparing coefficients yields 
Ci(w) = bi(w) for £ = 1, . . . , k. 
Let 

m 

q m (z,w) = S ^bi{w)z e for m G N 
and consider the Mellin transform of log <& w (x) — qk{x, w) 

oo 

(5.3) M w {s) = J (log$ w (z) -q k (x,w))x s ^dx. 

o 

The function M w (s) is holomorphic in the strip — k — 1 < < — p by general properties 
of the Mellin transform (cf. [7J. In order to obtain a meromorphic continuation of M w (s) 
we split the integral in (J5.3|) into two parts: 

l 

(5.4) AfW(s) = / (tog 0*0 -q k {x,w))x s - l dx 



o 

oo 



(5.5) M$\s) = Qog$ w (x)-q k (x,w))x'- 1 dx. 



i 

It is clear that Mjp(s) is holomorphic in 9fts > —k — 1. From the expansion (J5.2)) of 
\og$ w (x) we immediately get a meromorphic continuation of M^\s) into 3?s > — r for 
r > k + 1, r G N by the observation that 

i 



M w\ s ) = j {\og$ w {x) -q r -i(x,w))x s - 1 dx+ 







s + £' 

k+l 



The right hand side of this equation is a meromorphic function in > — r and therefore 
provides the meromorphic continuation of Mw\s). Thus Mw\s) has simple poles at 
s = —£ for £ G {k + 1, k + 2, . . .} with residues be(w). 

(2) 

We proceed similarly for Mi, (s). This function is holomorphic for < —p. From 
Theorem ^ we know that the periodic function F has a uniformly convergent Fourier series 



(5.6) F(u) = fn 



rneZ 
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From the fact that F{u) is holomorphic for \Qu\ < [3a (here and later we set a = we 
actually get that 

f m = O e {exp(-27ia((3 - e)\m\)) for all e > 0. 

Thus we can write 



00 / 

1 I log$ w (x) — q k (x, w) - fn 

■[ \ m£Z 



M%\s) = I [\og<S>M- qk (x,w)-Yj m x^ im ° + l ^ + log(^)x s - 1 dx 



\^ f m (\oga d A 1 

> — h logl— w) 



s + p + 2nima \ d — 1 

where the right hand side is a meromorphic function in C with simple poles at s = — p — 
2mmo~ (m G Z) with residue / m and at s = with residue — lo f_°^ — log(— w). This provides 

the meromorphic continuation of Mw\s). 
Thus we have proved 

Lemma 3. The function M w (s) given by (|5.3|) admits a meromorphic continuation to the 
whole complex plane with simple poles at s = —£ for £ G {/c + 1, k + 2, . . .} (with k = [p\), 
s = 0, and s = —p — 2-nima (m G %) with o = j^^- The residues at these poles are given 
by (the constants be(w) are given by ()5.2|) ) 

Res M w (s) = b e (w) for £ G {k + 1, k + 2, . . .} 

(5.7) Res M w (s) = - \og(-w) 

s=o d — 1 

Res M w (s) = f m for m G Z. 

Under the conditions of Theorem^ Corollary^ or Corollary \5\ infinitely many values f m 
are non-zero. 

Similar arguments applied to the function 

f log$(x) - logx - Yl beX J ^ dx 

with bi given by 

(5.9) log^ = f>* 

x i=i 

yield 

Lemma 4. The function Mq(s) given by ()5.8|) admits a meromorphic continuation to the 
whole complex plane with simple poles at s = —£ for £ G {k + 1, k + 2, . . .} (with k = [p\), 
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and s = —p — 2-nima (m G %) with a = and a double pole at s = 0. The residues at 
these poles (the principal part resp.) are given by (the constants hi are given by (|5.9|) ) 

Res M (s) = b e for£e{k + l,k + 2,...} 

s=—£ 

(s.io) Mo^--— - + •■■ 

Res M (s) = f m for m e Z. 

s=—p—2nimcr 

Under the conditions of Theorem^ Corollary^ or Corollary^ infinitely many values f m 
are non-zero. 

In order to use the information derived for M w (s) for the meromorphic continuation of 
the Dirichlet series 



(5.11) &,«,(*)= E » 



At>0 



we derive an alternative expression for M w (s) based on (|5.1jl 

oo 

M w (s) = / (\og$ w (x) - q k {x,w))x s ~ l dx 



(5.12) 



oo 



*«,(-m)=o o 

Furthermore, the second integral in ()5.12|) can be evaluated as 
(5.13) M w {s) = &,«,(-*)- " 



S Sin 7TS 



by the fact that 



oo 

2 







logd + .r) -ac+ — H h (~l) fc ^-l x'^dx = — - — for —k — l<^ts<—k 

2 k J s sin 7rs 



and the fact that the right hand side is the Mellin transform of a "harmonic sum" (cf. |2*H1 
29J). Equation (|5.1Hj) provides us with the meromorphic continuation of C$,ui( s )- This is 
somewhat a reversion of the ideas used in 
Summing up, we have obtained 

Theorem 7. Let $ be the solution of the functional equation (|3.1|) . Then the Dirichlet 
generating function C*,«;( s ) of the solutions of $(—//) = if for w < gwen m f!5.11|) 
admits a meromorphic continuation to the whole complex plane. Under the conditions of 
Theorem [3J Corollary ^| or Corollary [3] there exist infinitely many simple poles at points 
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of the form s = p + litimo (a — -^j, m G 1). The following special values are known 
(p = and f m are given by tfaty) 

Res C<s>,w(s) = — — (p + 2iximo~) sin rr (p + 2iximo~) for m G Z 
(<s>,w(rn) = for m < p,m G Z 

C*,«(o) = ^37 log (~ w ) 

C*,™( m ) = (-l) m_1 m6 m (w) /or m G N, 

where b m (w) is given by (|5.2|) . 

Similar arguments applied to §>(x)/x yield 

Theorem 8. Let $ 6e t/ie solution of the functional equation IjH.lj) . T/zen £/ie Dirichlet 
generating function C*,o(s) o/t/ie solutions o/$(— //) = /orp > gwen m (|5.11jl admits a 
meromorphic continuation to the whole complex plane. Under the conditions of Theorem^ 
Corollary ^| or Corollary [3] t/iere exzsi infinitely many simple poles at points of the form 
s = p + 2irimo~ (a = j^x> m E Z). The following special values are known (p = 



fm 

Res C*,o(s) — — — (p + 2irima) sin7r (p + 2 r nima) for m G Z 

C*,o(^) — /or m < p,m G Z \ {0} 
C*,o(0) = 1 

Ci,„(0)^ 

C*,o(m) = (-l) m - 1 m6 m /ormGN, 
where b m is given by (|5.9j) . 
Remark 3. By Lemma El we have 

OmO) = lim V (A n 

n— >oo ' » 

Zgp(-")({w}) 

which is the function ( PtW (2s) studied in jJTJH2|- 



z) \ 



6. Meromorphic continuation of the zeta-function 

Up to now we have only considered the eigenvalues originating from m = in (|2.1|) for 
fixed w G A. In general the multiplicity of eigenvalues for m > is given by a linear 
recurrent sequence (3 m {w). This sequence has a rational generating function 

oo 

B w {x) = ^ P m (w)x 
m=0 



P W (X) 
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with P w , Q w e Z[x}. Let r w denote the radius of convergence of B w , then by Pringsheim's 
theorem pole of B w ; let k w denote the order of this pole. Then the positivity 

of /3 m (w) implies that the coefficient Ck w {w) in the Laurent-expansion around 

1 - ^ 



is positive. 

Putting everything together, we can write 

(6.D aw= E ^E^P^w, 

—Au=/iu u>eA 

which provides the meromorphic continuation of £a to the whole complex plane. If log A d < 
\ds then all the functions (,<$>,w(s) are holomorphic in > \ds — £ for some s > 0. Since 
(a(s) has a simple pole at s = \d s by the fact that L(x) x x^ ds , at least one of the rational 
functions B w (x) has to have a pole at x — X~^ ds . All poles of functions B w at x = \~^ ds 
have to be simple. Let W denote the set of those w with r w = \~z ds . By our observation 
on the sign of C\{w) for w G W the Dirichlet series 

has positive coefficients. We apply [221 Theorem 9.5, p. 184] to the function i](s) to see 
that r](^ds + ikr) = cannot hold for fixed r and all k G Z \ {0}. Thus the function 

has a simple pole at s = ~ds and at least two non-real poles on the line 3?s = \ds- Since 
the remaining summands in (|6.1|) have no poles on this line, these points remain poles of 

CaOO- 

Summing up, we have proved: 

Theorem 9. Let G be a p. c. f. self-similar compact fractal, whose Laplace operator A 
admits spectral decimation in the sense of Definition^ with a polynomial of degree d. Then 
the Dirichlet generating function of the eigenvalues of A 



c*M= E T, 



has a meromorphic continuation to the whole complex plane with poles contained in a finite 
union of sets {pu + luima m 6 Z}, where a = and X is the parameter coming from 

spectral decimation. There is a simple pole at s = ^d$ (d$ denoting the spectral dimension 
of G). If dog A d < T^ds then Ca( s ) has at least two non-real poles on the line 3?s = \d$- 

Remark 4. The case of G = [0, 1] which gives the Riemann zeta function and has log A d = 
^ds shows that the condition log A <i < is needed for the last assertion. 
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Theorem 10. Let G be a p. c. f. self-similar compact fractal, whose Laplace operator A 
admits spectral decimation in the sense of Definition^ Then the following are equivalent: 

(1) Ca(s) has at least two non-real poles in the set ^d$ + I-kwHj (a = j^^j. 

(2) the limit\\m x ^ 00 x~^ ds L{x) does not exist, where L(x) denotes the eigenvalue count- 
ing function (jl.ljl 

(3) the limit lim^o+ P{t)t2 ds does not exist, where P(t) denotes the trace of the heat 
kernel (O). 



Proof. "((IJ)4^(j2J)": the general theory of Dirichlet series (cf. [TT]) implies that Ca(s)(A s — 
\z ds ) is of polynomial growth along vertical lines. Thus for some fceN the last integral 
in the Mellin-Perron formula (cf. |39|) 



E 



x) 2ni J s(s + 1) • • • (s + k) 



ids— e+ioo 



2 



E 



Res Ca(s)x s — 1 + / Ca(s)x s 



s =±d s +2-Kiia s(s + 1) • • • (s + k) 2ixi J s(s + 1) • • • (s + k) 

converges for some s > 0. The last integral is 0(x^ ds ~ £ ). The series constitutes x^ ds times 
a Fourier series in log A x, which is constant, if and only if there are no non-real poles on 
the line 3fts = \ds- Since the left-hand side of the equation is an iterated integral of L(x), 

x~^ ds L(x) shows oscillations, if and only if the right-hand side does. 
"(HJ)^©" : a g a i n by Mellin inversion we get 

lfifs+1+joo 



P{t) = in J ^(^(^ ds 



2<is+l — too 



2 



lds—e+ioo 



E 



2" a " 1 

Res ( A (s)T(s)t- s + / ( A (s)T(s)t- s ds. 

-d s +2mlcT 2m J 



^\d s +2-Kilo 

ij 

■^ds—e—ioo 

Again the series is t~z ds times a Fourier series in log A i and the last integral constitutes an 
error term 0{t-^ ds+£ ). The same argument as above proves the presence of oscillations, if 
and only if £a has non-real poles on the line 5is = \ds- □ 

Remark 5. In [T7|. p. 105] it is conjectured that the limit \im. x ^, oa x~^ ds L{x) does not exist 
for non-integer ds in the "lattice case" . The non-existence of the limit has been confirmed 
in the case of existence of localised eigenfunctions (cf. E3J 1221) an d i n the case of the 
existence of spectral gaps (cf. 
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Remark 6. In [H]] it has been proved that the n-step transition probabilities p n (x, y) (the 
discrete analogue of the heat-kernel) on certain self-similar graphs satisfies an asymptotic 
relation 

p n (x, y) ~ n~i ds F{\og x n) 

with a periodic function F. It is shown that F is not constant, if the Julia set of a related 
rational function is a Cantor set. In it is proved for the same class of graphs that F is 
not constant, if d$ is irrational. 

Remark 7. Theorems PI and [TU1 together confirm the conjecture (cf. fTJ p. 105]) that the 
limit lim^oo x~? ds L(x) does not exist for fractals admitting spectral decimation in the 
sense of Definition [T] with a polynomial of degree d and log A <i < \ds- 

7. Examples 

Example 1. In EH] the eigenvalues of the self-similar Dirichlet-Laplacian on the K- 
dimensional Sierpihski gasket have been identified. The i^-dimensional gasket admits 
spectral decimation with the polynomial Pk{%) — x(K + 3 + x). The eigenvalues of A are 
derived from the preiterates of the elements of the set A = {—2, —(K+ 1), — (K + 3)} with 
multiplicities 

1 for m = 1 
otherwise, 

(3 m {-{K + 1)) = ^^{K + l) m - 1 - for m > 2, 

P m (-(K + 3)) = ^y^(K + IT' 1 + for m > 1. 

This yields (with A = K + 3) 

oo 

(7.1) ( A (s) = A- s C*,- 2 (s) + C*.-(tf+i)00 E /W-(# + 1))*"™+ 

m=2 

oo 

U(A S -(K+1)) + 2(A--1) ) 

We first observe that the poles of C*,-( s ) at s = p + 2TTima, which are given by X s = 2, 
cancel by the fact that the residues of £§ }W at these points are independent of w and the 
values of the rational factors sum up to 0. Furthermore, we observe 

C#,-(*+3)(s) = - l)C*,o(s) 
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by the fact that $(V) = —{K + 3), if and only if &(\z) = and &(z) ^ 0. This implies 
that the poles at s = 27iima of the rational factor of the last summand in (J7.1)) are 
cancelled. Using Mellin transform, we see that the cancellation of poles on the line = 
is equivalent to the existence of the limit 



\og x (-w) oo / 1 \ 

hm - . n 1 - 



We did some numerical computations, which indicate that there are oscillations for w = —3 
and K = 2. Nevertheless, the pole of the rational factor at s = is cancelled by the zero 

Of (&-K-1- 

From this we derive, for instance, 

. . K + 1 / \og(K + 1) 
Ca(0) = — — [1-2- h 



Ci(0 ) . ,og 2 + 2 -^£^l(K + 1) - log{K + 3 ) 

+ K + 1 C" (0) - K + 1 C" (0) 

WD A ' 2 + 3A '" 1 



2(AT + 2)(AT + 3) 

Example 2. In jlu|HTlH2*] the spectrum of the Neumann-Laplacian on the 2-dimensional 
Sierpihski gasket has been studied. It has been shown that it admits spectral decimation 
with the polynomial p(x) = x(5 + x). The Neumann-eigenvalues are derived from the 
preiterates of the set A = {—3, —5} with multiplicities 

gm— 1 _|_ g 

An(-3) = for m > 1 

gm— 1 j 

/5 m (-5) = for m > 2. 

This gives 

oo oo 

Us) = (*M S ) E /U- 3 ) 5 ^ + C*.-s(s) E ^(- 5 ) 5_mS = 



m=l 



+ 77777" 77 ) C»,-3(a) + ( 77777^77 " 77777^7 ) S^-bOO. 



,2(5 s -3) 2(5* -1)/ ' w V 2 (5 s -3) 2(5 S - T 

which is in accordance with jU |32]- Notice that there is no pole at s — 0, the poles 
at the solutions of 5 s = 1 are cancelled in the second summand by the observation in 
Example ^ Numerical experiments seem to indicate that the poles of the first summand 
on the imaginary axis do not cancel with zeros of C$,-3- 

By the same arguments as in Example ^the poles at the solutions of 5 s = 2 are cancelled. 
Furthermore, in the second summand the poles at the solutions of 5 s = 1 are cancelled by 
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zeros of C*,-5- The poles of the rational factors at s = cancel and we obtain 

CA(0) = ^log 5 3-i 

>,, n x 3cg,- 3 (o)-q,- 5 (Q) , , 

C * (0) = 41oi5 l0g3 



Ca(1) 
Ca(2) 



7 

30 
1 

150 



The second derivatives ($ w {0) can be computed numerically by the following obser- 
vations. Around s = we can write M w (s) = — log (~ w > + H w (s) with a function H w 
holomorphic around 0, and by (j5.13|) we have 2H W (0) = C$ _ w (0). We have 

(1 - 2 • 5 s )H w (s) + rf lQg( ~ w) (5* - 1) = 

s 

!°g 7^ fxTTV^ ^ dx + / log X s - 1 dx, 



M 1 ~^^l5- 



and setting s = 



I oo 

\ dx f. ( 1 - z:$(a;) \ dx 



H w (0) = 2 hg(-w) log 5 - y log ( (1 _r $( , ))2 j - - y log [ 

o 5 1 



-^(l-^(f)) 2 



x 



The first integral can be computed as the rapidly convergent series 



log ( 1 ™ {X) ) ^(1 - 2 ■ 5-) 



The power series for log(l — —&(z)) has radius of convergence the modulus of the smallest 
solution of $>(z) = w, which is much larger than 1 in the cases of interest. This gives 
exponential convergence. 

For computing the second interval we observe that <&(x) > exp(0.9x p ) for x > 5, which 
can be shown by discussing $ on the finite interval 5 < x < 25 and the extending by 
the obvious inequality <3>(5x) > $(x) 2 . This together with an inequality for the logarithm 
yields (for T > 25) 

0> flogf )~^L 2 ) ->(2w-l) f°°eM-^)-. 
J T \-w(l-^(f))V x J T x 

Thus the improper integral can be computed by truncation at a finite value and estimating 
the remainder integral as above. 
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Using these ideas we computed (using Mathematica) 

F_ 3 (0) = 5.23995 51500... 
F_ 5 (0) = 9.06601 63789... 
Q0) = 0.96852 21499.... 

If there were no poles of Ca on the line 9ffe = 0, exp(— Ca(0)) were the value of the regularised 
product of the eigenvalues (cf.[12j), or the Fredholm determinant of A. 

Acknowledgment. The second author is indebted to Michel Lapidus and Alexander 
Teplyaev for helpful discussions. The authors are grateful to an anonymous referee for 
valuable suggestions. 
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